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We study a modified version of a U{1)b-l gauged MSSM tliat was recently sliown 
to produce a new source of leptogenesis through the CP asymmetry of sneutrinos and 
antisneutrinos [1]. By taking all superpotential terms and couplings between the MSSM 
Higgs and B — L scalar sectors into account we find that the model allows a large enough 
CP violation to explain the observed baryon number to entropy ratio. Monte Carlo 
analysis shows that a large amount of CP violation can be produced in the decays of the 
B — L Higgs bosons and that there are two dominating channels that drive CP violation. 

1 Introduction 

Baryogenesis via leptogenesis [2, 3j is one of the most appealing scenarios of explaining 
the observed excess of matter over antimatter in the Universe [4J indicated by the baryon 
number to entropy ratio 

— = (8.75 ± 0.23) X 10"^^ (1.1) 
s 

In these models right-handed singlet neutrinos, which are also responsible for the nonzero 
masses of Standard Model (SM) neutrinos j5], decay to SM leptons and the SM Higgs 
doublet creating CP asymmetry and violating lepton number. The resulting net lepton 
number is then converted to baryon number by sphalerons [6j. The basic picture of 
leptogenesis has been accommodated into supersymmetric (SUSY) models in various 
ways [3 m |9]. In [7], the MSSM superpotential is augmented with interactions between 
singlet right chiral neutrinos and MSSM lepton and Higgs superfields, and a Majorana 
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mass term is also included. With this setup sneutrinos and other superpartners can run 
in the loop diagrams and the decay products include the superpartners of SM leptons and 
Higgs bosons as well. In addition, sneutrinos can decay like their fermionic superpartners. 
References jHl [9] consider the effect of soft supersymmetry breaking to leptogenesis (soft 
leptogenesis) . It turns out that a single sneutrino generation can produce the required 
CP violation and net lepton number as opposed to the standard leptogenesis scenario. 

A natural way to extend SM by singlet neutrinos is to gauge B — L. Simple ex- 
tensions of SM with gauged B — L symmetry, i.e. containing the subgroup U{1)b-Li 
must accommodate three right-handed singlet neutrinos to cancel the triangle anomaly 
Gauging B — L within SUSY models would help understand R parity with 
the transformation R = (— 1)3(b~^')+2S' ^j^qj^ rjj^g breaking oi B — L can be attributed 
to additional Higgs fields that carry an even B — L charge and these Higgs bosons also 
generate the large Majorana masses for the right-handed neutrinos. 

A further modified model of supersymmetric leptogenesis was presented in [Ij where 
it was found that MSSM extended with gauged U{1)b-l gives rise to a new source of CP 
violation. Namely, the new heavy Higgs bosons that spontaneously break the gauged 
U{1)b-l symmetry can undergo decay into and A^* thus creating an asymmetry 
between these two. This new asymmetry is converted to conventional lepton asymmetry 
as the sneutrinos A^ decay into MSSM leptons and Higgs bosons and their superpartners. 
Both resonant leptogenesis (HJ-jl^ and soft leptogenesis arise in this model: CP violation 
is due to the complex parameters in the soft SUSY breaking sector and the heavy Higgs 
bosons are degenerate prior to the onset of soft SUSY breaking. After SUSY breaking, 
the Higgs boson masses receive suppressed contributions leading to a quasidegenerate 
neutral boson spectrum. 

We elaborate on the model presented in [1] by including in the superpotential and 
soft SUSY breaking potential the MSSM Higgs sector that couples to the field 5* and 
terms such as and AA. Thus, the sources of CP violation are the soft terms as 
well as the terms added to the superpotential, and as a consequence the model does not 
exhibit soft leptogenesis alone. The mass spectrum for the heavy Higgs sector is assumed 
to remain quasidegenerate and so resonant leptogenesis pertains to this model. In the 
present paper, we compute the CP violation parameter and lepton number with these 
modifications and investigate the allowed parameter regions with Monte Carlo methods. 
In section 2 we present the model, CP violation is studied in section 3 and the results of 
the numerical analysis are presented in section 4. The results are discussed in section 5. 
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2 The model 

The model is based on the gauge group SU{3)c x SU{2)lX U{1)y x U{1)b-l- The B-L 
charges for relevant chiral superfields A, A, 5, N, e'^, L, Q, u'^ and d!^ 
are +2, —2, 0, +1, +1, —1, +1/3, —1/3 and —1/3, respectively. The A^j fields receive 
Majorana masses through the vacuum expectation value (A) that breaks B—L symmetry. 
The most general superpotential for these fields reads then as 

1 

2' _ 
-Y^^'L^N.Hu + liHuHd + MiS^ + M2 AA + Y^S"" + Y^SHuHd, (2.1) 



W^-^ = XS{AA-M^) + -fijN,NjA 



the D term potential is 



VE-' = l{9' + 9"m'f + m'-\H',\^-\H-\r (2.2) 

+\9'\h:[ht + KH,r + 2^?i(iAp - \A\r 

and the soft SUSY breaking potential is 

Vf^j/' = [b{H+H^ -H^H°^) + aiS + hS^ + b2AA + ciS^ + C2SAA (2.3) 

+c,S{H:[H, - H^H',) + MlAN.N^] + h.c. 
+mljHu\'' + mj.jHdl'' + ml\S\'' + ml\A\'' + m|| Ap. 

The SU{2)l, U{1)y and U{1)b-l gauge couplings are denoted by g, g' and gB, respec- 
tively. For convenience, we move to unitary gauge by making the transformations 

A = -L(|M|+Ao)e''-^-^', 

A = ^(|M| + Ao)e-«^5^^'+''^^^ (2.4) 
v2 

Minimizing the scalar potential V = Vp~^ + V^'^ + V^^J^^ w.r.t. Re(i7„), Re(i7rf), 
\m.{Hu) and Im(ifrf), Re(S'), Im(S'), Re(Ao), Im(Ao) allows us to eliminate the parame- 
ters Tn^j^, '^'iidi = Re(— AM^), Re(ai), Im(ai), Re(62) and Li = Im(— AM^), respec- 
tively, and also, b is fixed. We are thus left with the free parameters Mi, M2, Yi, 1^3, 
bi, Ci, C2, C3, Af, m|, m\, and (S). The gauge coupling gB does not appear in any 
of the mass eigenvalues or mass eigenvectors so it is not included in the Monte Carlo 
analysis. 

Our minimization of the scalar potential departs from the corresponding procedure 
carried out in [T] in that we eliminate some of the soft parameters and randomize (S) in 
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the vicinity of 1 TeV. The vacuum expectawtion value (VEV) of S is nonzero, (S) 7^ 0, at 
SUSY breaking and the vacuum expectation value has to be of the SUSY breaking scale 
[T]. Because of the alterations made to the superpotential, (5*) is nonzero before SUSY 
breaking and altogether (S) ~ 1 TeV after soft SUSY breaking. The VEV (S) appears 
in the formulas of Li and Lj. as well as the soft parameters that have been eliminated in 
the minimization procedure. 

From the scalar potential V = Vp + Vd + Ko/t we determine the mass matrix in the 
basis (ReHu, ImHu, KeHd, ImHd, KeS, ImS, RcAq, ImAo). This computation is described 
in the Appendix. The expressions for the masses and eigenvectors containing the per- 
turbative corrections are very complicated and finding special limits is nontrivial. By 
limiting different parameters we would obtain different limits for the masses so there are 
no unequivocal expressions for the masses or the eigenvectors. The four quasidegenerate 
heavy B — L bosons of this system decay to sneutrinos and antisneutrinos generating CP 
violation, which we investigate in the next section. 



3 CP violation 

Since we are dealing with a quasidegenerate system in the heavy B — L Higgs bosons, 
we expect that the main contribution comes from the interference between the tree level 
decay diagram and mixing diagrams in Fig. [T] To compute the CP violation parameter 

TV / N , 



\ TV ' TV ^ 

(a) (b) (c) 

Figure 1: The relevant tree level and loop diagrams depicting the decay of the heavy 
mass states Xi to sneutrinos. Sneutrinos or their fermionic partners can run in the loop. 

we also need the corresponding diagrams that produce N* pairs. The diagrams of Fig. 
[T] arise from the scalar potential part 

^(3) = IN\'J2 F\fj\,X, + {NN J2 Fm^X^ + h.c.) (3.1) 

i i 

and the Yukawa interaction between the neutrinos and the heavy states Xj 

CNNX, = NNYp^X, + h.c., (3.2) 
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where the couphngs -FjvTVi given later. We have restricted ourselves to 

the case where there is only one neutrino generation A^. The sneutrino states and N* 
mix with each other and form the propagating mass eigenstates 

1 

1 



AT, 



Equation fl3.ip now becomes 



(3.3) 



(3.4) 



The definition for the CP violation parameter 

^ Ejr(X, ^ NN) - f{Xi ^ N*N*)] 
EJr(X, ^ NN) + T{X, ^ N*N*)] 

is employed and after the standard calculation we find 



Ml - Ml 



^ (Ml -MIY + n 



(3.5) 



(3.6) 
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with the absorptive part of the boson and fermion loops 

1 



327r 



167r 



'1- 



4m; 



[{MI-2MI)(yIYf + Y^Y, 



FJij 



The sneutrino squared masses and neutrino mass are 



Ml, 



M, 



V2 



\f\\M\. 



(3.7) 



(3. 
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The couplings are given by 




\f\'\M\n,,, 



(3.9) 



ffNi 




f {riji + insi) , 



X*e-'*\M\in5^-^ne^) 



F. 



F. 



+■ 



1 



i{e ^Fjyff- — e ^F^c.), 




X 



) 



The coefficients Uij relate the fields (ReiJ„, Iniif^, KeHd, ImHa, KeS, ImS, ReAo, IniAo) 
to the mass eigenstates Xj as shown in the Appendix. 

Our result for the fermionic part in the second line of the absorptive part of the loop 
(13. 7p is slightly different from that in [Ij. We have checked our result with the help of the 
optical theorem that suggests (13. 7p is correct. The correction —2M'^(YpYF + YpYp) is of 
the same order as —2Mjf(YpYpe~^''^f + Y^Y^e^*'^-'') so the modification is not significant 
in magnitude compared to the dominating term M'^_{YpYp + YpYp). 

In the couplings (13. 9 p the combination X{S) + M2 appears. This sum equals the VEV 
of S in [Ij if we minimize Vp w.r.t. S before SUSY breaking and keep the dominating 
terms ~ |M|2. The resuhing VEV of S is (S) = M2/A and after soft SUSY breaking {S) 
receives further corrections from the soft potential. This correction — (c2|Mpe~*'^/2 + 
a*)/(|Ap|Mp) to our (S) is the result they found in [Ij for (S) by keeping the largest 
terms ~ |Mp in the 5* derivative of the scalar potential. 

The source for CP violation does not solely lie in the soft SUSY breaking sector in 
our model because the complex couplings Fj^j^- include complex parameters from the 
superpotential as well. These new contributions come from the heavy B — L sector as 
well as the coupling part between the B — L and MSSM Higgs sectors, Y^SHuHd- Also, 
the soft sector is modified and, in particular, the scalar B — L and scalar MSSM Higgs 
sectors are again coupled via the term c^SHuHd- 
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Table 1: The parameters and their scanned values. 



Parameter 


bcannecl values 


1 /\/f 1 


U.i-iU.U lev 


IV 1 
l^ll 


U.Ui-iU 


W 1 
\J^3\ 


iU — iU 


\h 1 


u.ui-iu lev 


Cl,2 


n ni inn t^q\^ 
u.Ui-iU.u lev 


C3 


1 n~6 1 n— 3 TqX/' 

lu — lu lev 


1 A 1 


n 1 inn t^q\^ 
u.i-iu.u lev 


2 2 


n 1 inn ta^2 
u.i-iu.u lev 


|A| 


0.01-0.1 


\M\ 


10^ - 10^ TeV 


I/I 


10-5 - 10-3 


1(5)1 


0.1-10.0 TeV 


l/^l 


0.1-1.0 TeV 



4 Numerical results 

In determining the viable parameter regions we impose some conditions. First, the excess 
baryon number created by decays is given by [9, Ij 

— ~ -8.6 X 10-V, (4.1) 
s 

where the washout factor rj can be ~ 0.1 at most [H [161 El [E]. Thus, e ~ —10-^ is 
required. Second, the heavy Higgs decay rates must be smaller than the expansion rate of 
the Universe. Finally, we want to ensure that the system is quasidegenerate and demand 
that the heavy Higgs boson mass differences are a few orders of magnitude smaller than 
the heavy Higgs masses, Mx,, themselves. Also, negative mass squared values for the 
Higgs bosons and sneutrinos may arise and we have to filter these out as well on physical 
grounds. 

By performing Monte Carlo analysis on the system with the new superpotential pa- 
rameters and soft parameters in Table [1] 
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and the conditions 



- 10 



5 



< 



e < -10 
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(4.2) 



< 




< 



1000 



we obtain the plots (each with around 1800 points) in Figs. The scan ranges of the 
soft parameters shown in Table [T]set the SUSY breaking scale in the vicinity of 1 TeV and 
it may vary up to 10 TeV. The parameter ranges given in Table [1] all allow for sufficient 
CP violation and the corresponding plots thus show a uniform distribution of points. On 
the other hand, the plots with masses and CP violation show more structure. The third 
line in (14. 2p is an ad hoc assumption made to ensure the quasidegenerate nature of the 
heavy B — L Higgs boson masses. 

It is in fact the first condition on e of the set (14. 2 p that most heavily restricts the 
parameter sets. The heavy Higgs bosons have masses in the range ~ 10^ — 5 x 10^ TeV, 
Fig. |2j Fig. |2] also shows that in reality there are no points ruled out that would widen 
the gap between the large masses Mx^ beyond the third condition in (14. 2p . Without the 
CP violation parameter condition the allowed heavy Higgs masses would span uniformly 
the whole interval from 10^ to 10^ TeV. Also, the sneutrino and neutrino masses would 
span the interval of roughly 10~^ — 10'^ TeV without any restrictions. When the conditions 
(14. 2 p are imposed, the allowed range for sneutrino masses is ~ 0.1 — 100 TeV with some 
points also in the region above ~ 10^ TeV while the neutrino mass lies in the interval 
~ 1 — 100 TeV also with some points in the region ~ 1000 TeV, Fig. [3J The occurrence of 
points in the larger mass region is due to the behavior of e versus sneutrino and neutrino 
masses. The CP violation parameter |£:| becomes larger than the allowed interval in 
(14. 2 p when e.g. Mfj^ > 10'^ TeV. Once the sneutrino masses attain values ~ 10'^ TeV and 
higher, \e\ is reduced and some points pass the first condition of (14. 2p . Similar behavior 
is observed in the neutrino masses as closer to 1000 TeV some points appear. The Higgs 
masses m/^, tua and itlh sit near 130 GeV and 1 — 10 TeV, respectively, with itia < ttlh, 
Fig. m Especially the lightest Higgs boson has a mass in the vicinity of 130 GeV after 
(14. 2 p are imposed. Thus, our setup seems to favor the lower end of the estimated mass 
range of the lightest Higgs boson. (For a review of the Higgs boson masses in SM and 
MSSM, see [19j.) 

One remarkable difference to [Ij is the fact that |£:| tends to be too large instead 
of too small. In fact, \e\ tends to be close to unity in our model with heavy Higgs 
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masses ~ 10^ — 10® TeV, sneutrino masses ~ 10^ — 10^ TeV, neutrino mass ~ 10^ — 10^ 
TeV. Also there seems to be a lower limit to the scale of the CP violation parameter, 
kl ~ 10~^. The channels making the largest contribution to e are X2 — > A^A^ via X3 and 
X3 — )■ A^A^ via X2. This is due to the 1-2 orders of magnitude smaller mass difference 



\Mx2 ~ ^-Xal, Fig- 2(a), compared to the other differences which are essentially similar 
to Fig. 2(b), The reason why \Mx2 — ^Xal is a few orders of magnitude smaller than 
the other mass differences is not clear-cut due to the complexity of the mass formulas. 
This mass difference comes closest to the resonance conditions \Mx2 — ^Xsl ~ T [15j 
although this limit is not reached, see Fig. [5l 

To compare our setup with [IJ we performed Monte Carlo on the model in [IJ with our 
definition of e and found that the phase space is significantly more constrained without 
the MSSM Higgs sector and other interactions introduced into the superpotential. The 
yield of points that satisfy (14. 2p is the order of 10~^ or less in [Ij whereas the extended 
model has a success rate of the order of a few percent. In pj, \e\ takes values many orders 
of magnitude below 10~^ and attaining the correct magnitude becomes the threshold. 
Thus, including the MSSM Higgs sector and other allowed terms radically alters the 
characteristics of the B — L gauged MSSM. Some differences arise also in sneutrino mass 
Mff_, and consequently and M^r, Fig. [21 The 1-100 TeV range values are favored 
in the extended model in contrast to the ~ 1000 TeV and above values in [IJ. The 
sneutrino and neutrino masses are of the same order of magnitude as shown in Fig. [6l 
with differences determined by the soft parameters in Table [T]that set the SUSY breaking 
scale. 

The minimization procedure on the scalar potential leads to expressions for Lj and 
Lr and for the soft parameters Re(ai), Im(ai) and Re(62) where Im(y3) is located in 
the denominator and does not cancel. This l/Im(Y3) dependence survives to the mass 
matrix perturbation SM^^ (lA.ip and so the perturbative masses and eigenvectors also 
have ~ 1/Im(y3) dependence. By reducing Y3 all heavy Higgs mass differences become 
smaller but they do not approach the values given by the model in [1]. 

Moreover, we have checked that the use of perturbation theory in diagonalizing the 
8x8 Higgs boson squared mass matrix is legitimate because the perturbative corrections 
to the boson masses are relatively very small, namely \Mxi — TnQi\ < 10 TeV and for the 
lightest MSSM Higgs boson {Mx^ — rnorl < 0.01 TeV. These differences are obviously 
many orders of magnitude smaller than the masses themselves. 

The scans for Mi, M2, I3 and C3 have been restricted. We have limited and C3 to 



10 



very small values because these directly couple the MSSM and B — L Higgs sectors. This 
would mean that the MSSM Higgs masses, especially m^, could receive large ~ 1 TeV 
corrections. Expanding the range for Y3 and C3 up to ~ 1 and ~ 1 TeV, respectively, 
yields less points that satisfy (14. 2 p compared to the values in Table [T] so it is justified to 
stick to the lower end of the range. Allowing Mi and M2 to take values up to ~ M makes 
the use of perturbation theory no longer valid and the VEVs of A, A and S would be 
significantly altered. It is likely that {S) would no longer be ~ 1 TeV and instead would 
increase many orders of magnitude. This situation is not considered in our analysis. 

Overall, the phase space of the model is rather complex and no single parameter from 
the superpotential or soft potential dominates the behavior of physical quantities. Thus, 
obtaining clear relations between physical quantities, e.g. \Mx2 — Mx^^l versus Mx4, is 
possible only if the parameters in the superpotential and soft potential are fixed to values 
that satisfy (14. 2p . Some parameters, e.g. Y3, C3, Yi and Ci, can be allowed to fluctuate 
around the values that pass ( 14. 2 p but this causes the masses Mx^ themselves to vary by 
very little, (9(10^^) TeV. Likewise, limiting physical quantities, e.g. Mj^r, instead of the 
parameters does not constrain the phase space enough to produce meaningful relations 
depicting the dependence between other physical quantities like masses, mass differences 
or the CP violation parameter. On the other hand, the Monte Carlo analysis unveils 
well the behavior and allowed regions of the MSSM Higgs boson masses corrected by the 
B~L sector, sneutrino masses, neutrino mass and CP violation parameter. Especially it 
has revealed that the lightest Higgs boson mass rrih seems to favor the ~ 130 GeV region, 
and that we have two clearly dominant channels and that these channels can produce a 
large amount of CP violation e < 1/2. 

5 Discussion 

We have found that matter-antimatter asymmetry via an asymmetry between and A^* 
can be produced in the B—L gauged MSSM with reasonable parameter values. We do not 
have to resort to finetuning because including the coupling term between the B — L and 
MSSM Higgs sectors and other allowed terms in the superpotential allows the model to 
produce a large excess of matter over antimatter. The situation is opposite in [Ij where \e\ 
tends to be many orders of magnitude below the required lower limit. The success of our 
model is largely due to the fact that two of the heavy Higgs masses and Mx-^ come 
very close to each other compared to the other masses and subsequently the resonant 
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condition becomes closer for these two particles. The decay channels X2 — > — > A^A^ 
and X3 — 7- X2 —7- NN thus dominate the CP violation parameter which can become 

kl < 1/2. 

The results from our model suggest that including the MSSM Higgs sector as well 
as other allowed interactions may change radically the characteristics of the physical 
system, in this case the amount of CP violation. Also, many other models that couple 
the MSSM Higgs bosons to the new Higgs bosons that spontaneously break the symmetry 
introduced beyond MSSM, see e.g. [201 EH l22], could be worthwhile studying in terms 
of leptogenesis. 

Acknowledgments The work of H.K. was supported by Finnish Academy of Science 
and Letters (Vaisala fund). 



A Diagonalization of the mass squared matrix 

In solving the neutral Higgs boson mass spectrum by diagonalizing an 8 x 8 matrix we 
rely on perturbation theory that is familiar from quantum mechanics (see e.g. |23|). 
The zeroth-order squared mass matrix consists of the MSSM Higgs block situated in the 
upper left corner and the B — L Higgs sector in the lower right corner. Modifications 
compared to [IJ such as the superpotential coupling Yi are included in the unperturbed 
B — L block. The corrections to the squared masses and eigenvectors are computed using 
the standard formulas 

^V„ = V<iS«kipi, (A.1) 



where Vn denote the eigenvectors and the eigenvalues of the zeroth order squared 
mass matrix and M'^j^q is given later. The perturbation to the zeroth order squared mass 
matrix is denoted by SM'^^^ and it equals the difference between the squared mass matrix, 
which is derived by differentiation w.r.t. Higgs scalar fields from the scalar potential V, 
and (lA.Sp . The perturbation thus holds various terms including the added couplings to 
the superpotential Mi, M2, Yi and 1^3 and their soft counterparts. We find the mass 
squared matrix for the scalar bosons in the basis 

XI = (Reif^t, ImiYu, Reifd, Imifrf, ReS", ImS", ReAo, ImAo) and obtain the following de- 
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compositions in terms of the mass eigenstates 



(A.2) 



where the coefficients riij are normahzed lengthy expressions including first order per- 
turbative corrections. The masses of the Higgs boson sector are denoted by and 
Mx2, are the heavy states that arise from the S — L sector. The hght Higgs states 

corresponding to the MSSM Higgses are denoted by = A, Xt = h and Xg = H. Then 
the Goldstone boson is massless Mx-^ = 0, the lightest Higgs h has mass Mx^, Mxg cor- 
responds to the mass of A and the H boson has mass Mxg- 

The zeroth-order mass squared matrix is in the basis (Rei^u, Inii^u) ^^H^, ImHd, 
ReS, ImS, ReAo, ImAo) 



M^BLO 



Vu 



-b- 
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bvd 

Vu 



b 








^ 2 



2 "T ». . 











b 



bvu 













as 








hs 



bs cs 














o\ 







bn 



(A.3) 



bn CD J 



QUYu + (Rc(|M| + (Ao))' + Im(|M| + (Ao))')(y2' + Yl) 

+3Fi,(Re(|M| + {^o) fY2^ - Im(|M| + (Ao))'F2i 

+2Re(|M| + (Ao))Im(|M| + (Ao))r2r-) cos(0) 

+3yH(-2Re(|M| + (Ao))Im(|M| + (Ao))y2i 

+Re(|M| + (Ao)) - Im(|M| + (Aq)) V2,) sin((/.) 

3ri,(2L, + (Re(|M| + (Ao))'r2. - Im(|M| + (Ao))'y2^ 

+2Re(|M| + (Ao))Im(|M| + (Ao))F2.) cos((/.) 

+(-2Re(|M| + (Ao))Im(|M| + (Ao))y2i 

+Re(|M| + (Ao)) V2, - Im(|M| + (Aq)) V2,) sm((/.)) 

-QL.Yu + (Re(|M| + (Ao))^ + Im(|M| + {l^o)f)iY^. + Y^) 

-3rH(Re(|M| + (Ao))'r2^ - Im(|M| + (Ao))V2, 

+2Re(|M| + (Ao))Im(|M| + (Ao))r2r) cos(</.) 

+3yH(2Re(|M| + (Ao))Im(|M| + (Ao))y2i 

-Re(|M| + (Ao)) V2, + Im(|M| + (Aq)) V2,) sin(0) 

^((5Re(|M| + (Ao))^ + 3Im(|M| + {^,)f){Yi, + Y^ 

+4L,F2.cos(0) + (Im(|Af I + (Ao))(-r2. + Y^r) 

+Re(|M| + (Ao))(r2i + r2,))(Re(|M| + (Ao))(F2i - Y^^) 

+Im(|M| + (Ao))(y2i + Y2r)) cos(20) + AL^Y^r sm(0) 

+2(-Im(|M| + (Ao))y2i + Re(|M| + (Ao))y2r-)(Re(|M| + (Ao))y2i 

+Im(|M| + (Ao))F2r)sin(20)) 

^(y2rcos(0) - r2i sin(0))(2L, + (Re(|M| + {^^)fY2i - Im(|M| + (Ao))^ 

+2Re(|M| + (Ao))Im(|M| + (Ao))y2r-) cos((/)) 

+(-2Re(|M| + (Ao))Im(|M| + (Ao))y2i 

+Re(|M| + (Ao))2r2r - Im(|M| + (Ao))'r2r) sin(0)) 

l((3Re(|M| + (Ao))^ + 5Im(|M| + {/\,)mY^, + Y^ 

-AUY2iCos{(t>) - (Im(|M| + {^o)){-Y2i + Y^r) 

+Re(|M| + (Ao))(y2^ + l2.))(Re(|M| + (Ao))(y2i - I2.) 

+Im(|M| + (Ao))(y2i + l2r)) cos(20) - 4L,y2r- sin((/.) 

+2(Im(|M| + (Ao))F2. - Re(|M| + (Ao))F2r)(Re(|M| + {Ao))Y2i 

+Im(|M| + (Ao))F2r)sin(20)). 
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In the above formualas, Y2 = A. The subscripts i and r denote the imaginary and real 
parts, respectively. The eigenvectors of (1A.3P are 

V, = J: ^ (0,-^;„,0,^;,,0,0,0,0)^ 

1 



2 



X (0, 0, 0, 0, 0, 0, ao - ^JAbl + [an - cd? - cd, 26^) 
1 



{ao + ^^hl + {ao - cdY - cdY + 46^ 



D 



X (0, 0, 0, 0, 0, 0, an + ^J ^hl + (a^ - cdY - cd, 2bD] 
1 

{as - v/46| + (as - csY - csf + 46| 



X (0, 0, 0, 0, as - \lAhl + {as - csf - cs, 2bs, 0, 0)^ (A.4) 
1 



{as + V46| + {as - csY - csY + 46| 



X (0, 0, 0, 0, as + \Ubl + {as - cs)^ - cs, 2bs, 0, 0)^ 



T 



= ^ (0,i;,,0,t;.,0,0,0,0) 

V7 = (cos(«),0,-sin(a), 0,0, 0,0,0)^ 
1^8 = (sin(a),0,cos(a), 0,0, 0,0,0)^ 
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with 



-{-8big" + g')vdvM - <f + (2& 

+(^'2 + g'')vdv^f{vl + vlffl-')l{{{g'^ + g'')vd{vd - v^)v^{vd + v^) 
-r1h{-vl + vl) - {-8b{g'' + - + (26 

+{2vdVu{2b + {g'^ + g^)vdVu)ff'^)/{{2vdVu{2b 

+2bi~v', + vl) - i-Sbig" + g')vdVu{vj - vD' + {2b 

H'^VdVu{2b + {g'' + g')vdv^))^)'/')) (A.5) 



and the squared masses are 




(A.6) 



{2b{vl + vl) + {g'^ + g^)vdVu{vl + vl) 



-^-8b{g'^ + g^)vdVu{vj - vl)^ + (26 + {g'^ + g^)vdVu)'{vj + vl)^) 
-^{2b{vl + vl) + {g'^ + g^)vdvM + vl) 



+^-8b{g'^ + g^)vdVy{vl - vl)^ + {2b + {g'^ + g')vdVy)^{vj + vl)^). 
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Figure 2: The mass difference \Mx2 — Mx^l versus in Fig. 2(a) and |M 



M 



versus Mx^ in Fig. 2(b) The difference |Mx2 — | is a few orders of magnitude smaller 
than |Mx4 — Mxr,\ and other mass differences. \Mx2 ~ ^x-i\ occurs in the amplitude of 
X2 — ?► X3 — 7- iViV which then dominates the CP violation parameter. 
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Figure 3: CP violation parameter versus neutrino mass Mjv- Both cases [T] (crosses) 
and the modified model (circles) are shown. The CP violation parameter \e\ tends to the 
higher end compared to \X\ ■ The sneutrino masses M^y^ spread around the Mn values by 
the sphtting due to soft SUSY breaking (13.81) satisfying the hierarchy < M^r < . 
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Figure 4: The CP violation parameter vs m/i in Fig. 4(a) and niH in Fig. 4(b) The 
masses rrih and m/^ contain corrections from the B — L sector on top of the predictions 
of MSSM. 
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Figure 5: The mass difference \Mx2—Mx:i \ versus the total decay rate of particle X3. This 
plot illustrates how the situation relates to the resonance condition \Mx2 — ^Xal ~ T- 




Figure 6: Sneutrino mass vs neutrino mass M^. The mass differences are of the 

SUSY breaking scale set by the soft parameters in Table [H 



